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Осуществлен синтез робастных алгоритмов управления нестационарными электромеханическими объектами, произведена 
оценка их чувствительности. 


Введение 

Принципиально важным свойством адаптивных 
методов обработки информации является робаст- 
ность, под которой понимается статистическая на- 
дежность метода и нечувствительность результатов к 
изменениям условий наблюдения. Основной осо- 
бенностью рассматриваемых объектов является их 
нестационарность: параметры объектов с течением 
времени изменяются в определённом диапазоне. Не- 
стационарность вносит принципиальные трудности 
как в изучение структурных свойств объекта (устой- 
чивости, управляемости и наблюдаемости), так и в 
разработку алгоритмов оценивания и управления. 

1. Описание модели объекта управления 

Пусть математическая модель функционирова- 
ния электромеханического объекта описана систе- 
мой линейных нестационарных дифференциаль- 
ных уравнений вида [1]: 

х(() = А(і)х(і)+Ь(()и(і), х(1 0 )=х 0 , (1) 

где х(1 ) - я-мерный вектор, компоненты которого 
определяют состояние электромеханического 
объекта в момент времени /; и(?)_- т-мерный век- 
тор управляющих воздействий; А(і) - матрица па- 
раметров объектов управления размерности (яхя), 
Ь(Г) - матрица влияния управляющих воздействий 
размерности (яхт); х 0 - начальное состояние элек- 
тромеханического объекта в момент времени / 0 . 

Так как сигнал с ПЭВМ поступает в дискретные 
моменты времени, дискретная модель, соответ- 
ствующая непрерывной модели ( 1 ), имеет вид: 

х(к + 1) = А(к)х(к)+Ь(к)и(к), х(0) = х 0 , (2) 


где матрицы А(к ) и Ъ(к) рассчитываются по форму- 
лам [2]: 




АГ А(і к У 


т=0 
оо А +т 


т\ 




( 3 ) 


Если обозначить через А(к) и Ь(к) суммы ./V пер- 
вых членов рядов (3) соответственно, то матрицы 
Д ѵ (4) и Ь^к) аппроксимируют А(к) и Ъ(к) с погреш- 
ностью_порядка 0(А/ ЛГ ). При этом число слагаемых в 
А(к) и Ъ(к ) можно задавать заранее (вычисление с 
фиксированной точностью) или определять авто- 
матически с помощью соотношения: 
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где е выбирается из условия обеспечения макси- 
мальной точности вычислений на ПЭВМ конкрет- 
ного типа. 

Такт к в (2-4) соответствует моменту времени 
4 =( 0 +кА1, А I- период дискретности по времени. 

2. Структура широтно-импульсного модулятора (ШИМ) 

В качестве управляющего воздействия рассма- 
тривается выходное напряжение импульсного пре- 
образователя, модулированного по широтно-им- 
пульсному закону и поступающее на исполнитель- 
ный элемент электромеханического объекта. Ши- 
ротно-импульсный модулятор (ШИМ) описан сле- 
дующими уравнениями: 
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и(к + 1) = 


'К й (к) при й (к) < — , 


при й (к) > 


9? 


Щк) = 


V п при 1 0 +кТ <к<1 0 + кТ+и(к), 

0 при 1 0 +кТ + и(к) < к < ( 0 + (к + \)Т, 


где и(к) - входной сигнал в такт времени / 0 +М/, 
ІН и Т- коэффициент передачи и период импульс- 
ного преобразователя электромеханического моду- 
ля соответственно, ІІ П - напряжение питания ис- 
полнительной части объекта управления. 


3. Синтез регуляторов функционирования 
системы управления объектом 

Для увеличения качества процесса управления 
объектом (увеличения быстродействия и точности 
достижения заданного значения) синтезируем нес- 
колько типов регуляторов. 

3.1. Пропорциональный регулятор 
на основе функций Ляпунова 2-го порядка 

Параметры регулятора будем определять, исхо- 
дя из минимизации квадратичного интегрального 
критерия: 

/ = | е т 0) (2 е(1) Рі , 

1 0 

где е(1)=у, ад (1)-у(1) - ошибка регулирования, 
т - знак транспонирования, 0 - диагональная ма- 
трица (яхя), элементы которой определяются ана- 
логично с п. 3.2. 

Для получения необходимой плавности движе- 
ния в окрестности заданной точки позиционирова- 
ния и асимптотической устойчивости возмущён- 
ной системы управления синтез управления при- 
водом может быть осуществлён на основе второго 
метода Ляпунова. 

Используя дискретное представление объекта 
(2) и функцию Ляпунова квадратичной формы: 

Ѵ(х(к)) = х(к) Г Рх(к), (5) 

где Р - положительно-определённая (яхя) - ма- 
трица, управление и(к) будем находить из условия 
минимума первой разности функции Ѵ{х{к))\ 

А Ѵ(х(к)) = Ѵ(х(к + \))-Ѵ(х(к)), 

вычисленной на решении системы (2). 
Минимизируя (5) по и(к) получим: 
й(к) = -((Г(к) т РВ(к)) т + 
+В(к) т РР(к))(В{к) т РВ{к)У\ 

где Р(к) и В(к) - линеаризованные по методу Эйле- 
ра в первом приближении матрицы: 

Р(к) = I + А(к)Т, В(к) = Р(к)Ь(к). 


Матрица Р определяется из условия минимума 
критерия (6): 

Р(к) т РР(к)-Р = -(). 

Построенный закон управления обеспечивает 
наиболее быстрое убывание функции Ляпунова 
(5), определённой на решениях системы (2). 

3.2. Пропорциональный регулятор 

на основе функций Веллмана 

Пусть модель объекта задана в виде (2). Будем 
считать, что вектор состояния х измеряем полностью. 

Функционал качества имеет вид [3]: 

Р = ^(х(к) т д{к)х(к) + 

к-0 

+и(к) т К и (к)) + х(2Ѵ) г Л х(2Ѵ), (6) 

где Л - неотрицательно определенные симметри- 
ческие матрицы размера (яхя), К - положительно 
определенная симметрическая матрица Тре- 

буется найти управление и(к) с полной обратной 
связью, минимизирующее функционал (13). Ура- 
внение Веллмана принимают вид: 

х(к) т Ох(к) + й(к) т К(к)й(к)+ ^ 

+В(к + 1, А(к) х + В(к) й (к)) 

Функция Веллмана ](к, х) определяется по фор- 
муле: 

Лк,х) = х т (к)Р(к)х(к), (8) 

где Р(к) - неизвестная неотрицательно определен- 
ная симметрическая матрица размера (яхи). Полу- 
чаем, что в задаче (6-8) оптимальное управление 
определяется соотношением: 

й (к) = К(к)х(к), 
к = 0,1,.. „(М-1), 

где К(к) - матрица коэффициентов усиления регу- 
лятора размера (дхп): 

К (к) = [К + В(к) т Р(к + 1) В(к)Т' х 

хВ(к) г Р(к + 1)А(к), к = 0, 1, (2Ѵ — 1), 

а матрица Р(к) размера (яхя) удовлетворяет уравне- 
нию: 

Р(к) = () + К(к) т КК(к) + 

+(А(к ) - В(к) К(к)) т Р(к + 1) (А(к) - В(к) К (к)), 

к = {Ы- 1) 1,0, 

Р(А 0 = Л. 

Матрица 0 оценивает отклонения реальных пе- 
ременных состояния от желаемых и представляет 
собой своеобразный “штраф” за большие ошибки 
при любом /„</<?;,. Матрица К, будучи всегда поло- 
жительной, оценивает стоимость управления. Фи- 
зически она характеризует затрагиваемую энергию 
на управление. Матрица Л характеризует ошибку 
управления в конечный момент времени. 


У(А„х) = тіп 
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Технические науки 


Предварительный выбор значений элементов 
весовых матриц 0 и К может быть осуществлен с 
помощью следующих рекомендаций: 

- обычно матрицы 0 и К назначаются постоян- 
ными и диагональными, т.е. матрица 0 содер- 
жит п ненулевых элементов ^ и , і=\,2,...,п, а ма- 
трица К— г содержит г элементов 

— принимаем, что максимально допустимые от- 
клонения переменных состояния х(к) в любой 
момент времени вносят в функционал качества 
одинаковый вклад. 

Распространяя аналогичные рассуждения и на 
отклонения сигналов управления и(к), можно 
записать: 


Ча = 




ѵ 1тах 


Яп, г п = 




и 1 тах 
1^7 тах ^ 


41- 


( 9 ) 


Здесь х ітах и Х/ тах — максимально допустимые от- 
клонение, определяемое техническим заданием, 
/-ой переменной состояния (/=2,3,..., п) и /-ого 
сигнала управления (/'=2,3,..., г) соответственно. 
Общий вклад максимально допустимых отклоне- 
ний переменных состояния должен приблизитель- 
но соответствовать общему вкладу максимально 
допустимых отклонений сигналов управления. 

У, Чіі Х ‘ тах = У Г л и У тах \ ( 10 ) 

1=1 і = 1 


- произвольно выбираем значение г п (например, 
г п =1) и по формулам (9, 10) вычисляем значе- 
ния остальных коэффициентов. 

Полученные значения весовых коэффициентов 
следует рассматривать как начальные оценки. Если 
отклонения по всем переменным достигают своего 
максимума не одновременно, формулы (9, 10) не- 
верно отражают требования, предъявляемые к си- 
стеме. Поэтому окончательный выбор весовых ко- 
эффициентов целесообразно производить после 
нескольких пробных процедур синтеза и модели- 
рования системы управления. 


3.3. Пропорционально-интегрально-дифференциальный 
регулятор (ПИД) 


Наиболее распространенными алгоритмами яв- 
ляются пропорционально-интегральный (ПИ) и 
ПИД алгоритмы цифрового управления. При пра- 
вильной настройке эти алгоритмы обеспечивают 
достаточно хорошее качество управления для боль- 
шинства объектов промышленной технологии. 


Рассмотрим процедуру вывода алгоритма ци- 
фрового ПИД-регулятора из соответствующего не- 
прерывного закона, имеющего вид [4]: 


й (/) = К 


е(/) + уг|е(/)Л 


+ Т, 


с2е(() 

Л 


(П) 


где К р , Т„ Т л - пропорциональный, интегральный, 
дифференциальный коэффициенты регулятора со- 
ответственно, е(і) - ошибка рассогласования. 


Запишем ( 1 1 ) в конечных разностях путем заме- 
ны: (=кТ ь А/ - период дискретности во времени: 


й(к) = К р 


Ф) + чг^Ф) г к +Т Л 


е(к) -е(к -\) 


т к 


где кх= 1,2,3,...,(/Ѵ-И), Т к - номер и величина перио- 
да квантования. 

Отметим, что при достаточно малых периодах 
квантования цифровой ПИД закон управления 
обеспечивает почти такое же качество процессов 
управления, что и исходный непрерывный закон 
(11). Для того, чтобы эффект квантования по вре- 
мени мало сказывался на динамику системы ци- 
фрового регулирования, рекомендуется выбирать 
период квантования из соотношения [5]: 


Т Т 

—А < Т <^~ 
15 * 5 ’ 


где Т 95 - время достижения выходным сигналом 
уровня 95 % от установившегося значения при по- 
даче на вход объекта сигнала ступенчатой формы. 

Настройку ПИД-регулятора (нахождения пара- 
метров К р , Т и , 7Д, будем осуществлять, руковод- 
ствуясь заданной относительной степенью затуха- 
ния (р [5]. Определяем значение логарифмического 
декремента затухания т: 

— 1п(1 -ср) 
т = — — . 

2 к 

Система автоматического управления будет 
обладать требуемой относительной степенью зату- 
хания, если амплитудно-фазовая характеристика 
разомкнутой системы автоматического управления 
будет проходить через точку на комплексной пло- 
скости (—1,70) [6]: ^=0,409, 7^=8,696, 7>=0,035. 


4. Оценка робастности 

Проведём анализ синтезированных регуляторов 
на робастность (реакцию на изменение параметров 
объекта). Будем считать, что объект является неста- 
ционарным, т.е. его параметры изменяются во вре- 
мени. Положим, что объектом управления является 
двигатель постоянного тока ДК-1. Переходные про- 
цессы по току двигателя без вынужденных возмуще- 
ний и при увеличении момента сопротивления в мо- 
мент времени 1 к =1 9 +к-\1, 7=200 (заданное значение 
по скорости равно 1) изображены на рисунке. 


Таблица. Характеристики синтезированных пропорцио- 
нальных регуляторов 


Тип регулятора 

Перерегулиро- 
вание, % 

Установивше- 
еся значение 

Время, такты 

Без воз- 
муще- 
ний 

С на- 
груз- 
кой 

Без воз- 
муще- 
ний 

С на- 
груз- 
кой 

Без воз- 
муще- 
ний 

С на- 
груз- 
кой 

К=1 

1,532 

- 

0,999 

- 

144 

- 

По функциям 
Ляпунова 

0 

0 

0,995 

0,983 

126 

300 

ПИД 

0 

0 

0,994 

0,987 

140 

261 

По функциям 
Веллмана 

0 

0 

0,995 

0,984 

122 

295 
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По функциям Ляпунова 

пид 

По функциям Б елмана 


Рисунок. Переходные процессы по току на обмотке двигателя, обеспечиваемые синтезированными регуляторами ( без вынуж- 
денных возмущений - слева и при увеличении нагрузки на валу в 3 раза - справа) 


По результатам, приведённым в таблице, можно 
заключить, что робастность регуляторов, синтези- 
рованных по предложенным алгоритмам, опреде- 
ляется в соответствии с конкретно поставленными 
задачами. В зависимости от требований, предъя- 
вляемых к объекту управления, поставленными пе- 
ред ним целями и возможными затратами, может 
быть выбран тот или иной алгоритм управления. 

В нашем случае наиболее подходящим является 
ПИД-регулятор. 
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Заключение 

В современной теории управления нестацио- 
нарными электромеханическими объектами боль- 
шое внимание уделяется робастности алгоритмов к 
действующим возмущениям. В данной статье были 
исследованы на нечувствительность несколько ти- 
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